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We investigate irreversible aggregation in which monomer- monomer, monomer-cluster, and 
cluster-cluster reactions occur with constant but distinct rates Kmm, Kmc, and Kcc, respectively. 
The dynamics crucially depends on the ratio 7 = Kcc /Kmc and secondarily on e = Kmm /Kmc- 
For e — and 7 < 2, there is conventional scaling in the long-time limit, with a single mass scale that 
grows linearly in time. For 7 > 2, there is unusual behavior in which the concentration of clusters 
of mass k, Ck decays as a stretched exponential in time within a boundary layer k < k* oc t 1 ' 2 ^ 1 
(k* oc hit for 7 = 2), while Ck oc t~ 2 in the bulk region k > k* . When e > 0, analogous behaviors 
emerge for 7 < 2 and 7 > 2. 

PACS numbers: 02.50.-r 



I. INTRODUCTION 

In this work, we investigate a simple aggregation pro- 
cess in which three kinds of reactions occur with constant 
but distinct rates: 





Mi 


A 3 , 




A x - 




1 . 
— A j+1 


3>2, 


Ai- 


lA 3 


1 , A.,, 


i,3>2 



Equivalently, the reaction matrix K+j has the value e in 
the upper left corner, the value 1 for all elements along 
the top and left edges of the matrix, and the value 7 in 
the rest of the matrix. 

Despite its simplicity, this model exhibits rich dynam- 
ics in which both scaling and universality can be violated. 
The model was initially studied by Hendriks and Ernst 
PJ to describe polymerization in which addition processes 
(reactions involving monomers) occur more readily than 
aggregation (reactions between i and j, with i,j > 2). 
They found evidence of unusual kinetic behavior by an 
implicit solution for the cluster concentrations; see also 
for related work. Our focus is on the opposite limit in 
which aggregation dominates over addition. By explicit 
solutions of the rate equations, we elucidate the full range 
of behaviors for this system. 

From a general perspective, the origin of the anoma- 
lous behavior in addition/aggregation stems from the fact 
that the simplest case of aggregation with a constant re- 
action rate is actually a marginal system. As discussed 
by van Dongen and Ernst Q , aggregation can be broadly 
categorized by whether reactions among large clusters 
dominate (type I systems, in their nomenclature) or 
whether reactions between large and small clusters dom- 
inate (type III systems). The former generically leads to 
monotonically decaying cluster mass distributions, while 
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the latter gives peaked distributions. The constant reac- 
tion rate system (type II system) is marginal by being at 
the boundary between these two behaviors. 

Because of this marginality, the kinetics of constant- 
kernel aggregation is sensitive to small perturbations 
in the reaction rates. Previously-studied examples of 
this feature include aggregation involving two distinct 
monomeric units, so that the reaction rates between 
two even masses, two odd masses, or an even and an 
odd mass are naturally distinct. Here non- universal and 
non-scaling behavior arises as a function of these rates 
0, 0. Another example is aggregation with the rate 
Kij = 2 - q i - q j with < 9 < 1 |. Although this 
reduces to a constant- kernel system for i,j — > 00, the 
g-dependent terms lead to unusual kinetics. The model 
studied here is in the spirit of the first example, except 
that we perturb only an infinitesimal fraction of the re- 
action rates in a constant-kernel system. It is striking 
that such a small change in the reaction rates has such a 
profound influence on the kinetics. 

In the next two sections, we focus on the special case 
of "sterile" monomers, where monomer-monomer reac- 
tions do not occur. In Sec. II, we first determine the 
monomer and dimer concentrations and show that dif- 
ferent behaviors arise for 7 > 2 and 7 < 2. Then in 
Sec. Ill, we derive asymptotic results for the cluster con- 
centrations. We then study, in Sec. IV, these same three 
cases when monomer-monomer reactions can also occur. 
Sec. V gives a summary as well as a discussion of the 
equivalence between the model and a diffusion-controlled 
process in which monomers have different diffusivity than 
other clusters. 



II. STERILE MONOMERS 

When monomers do not interact among themselves 
(Kn = 0), we shall show that the system exhibits three 
distinct kinetic regimes, with unusual time dependences 
and breakdown of conventional scaling in two of these 
cases. 
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The rate equation for the density Ck(t) is 7] 

«*(*) = J E K >Mt)cAt) - Cfc(t) E *kj c j(*). ( 2 ) 



i+j=k 



J>1 



where the overdot denotes time derivative. The first 
(gain) term accounts for the formation of clusters of mass 
fc (fc-mers) as a result of the bimolecular aggregation of 
i-mers with j-mers, with i + j = k. The second (loss) 
term accounts for the aggregation of the fc-mers with 
any j-mer, thus leading to a decrease in the fc-mer con- 
centration. For the processes given in and with the 
additional constraint of Ku = e = 0, the explicit rate 
equations are 
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-ciR, 

-c 2 [ci + jR] 
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where we drop the time argument for compactness and 
introduce the shorthand R = X)fc>2 c k for the concentra- 
tion of reactive clusters (those with mass > 2). Notice 
also that there is no production of monomers and dimers 
in the reaction; they can only disappear. 

The governing equation for R is obtained by summing 
Eqs. © to give R — — ^R 2 , whose solution is 

R = R{0)t-\ T =l + -jR(0)t. (4) 

With the explicit expression for R, we integrate the equa- 
tion for monomers to give 



Cl = Cl (0)r- 2 ^. 



(5) 



The crucial point is that the value of 7 controls which 
of R or ci decays more rapidly. For 7 < 2, monomers 
are asymptotically irrelevant and conventional constant- 
kernel aggregation kinetics arises. For 7 > 2, however, 
the large disparity in the reaction rates K\j and K%j leads 
to a long-lived residue of monomers and these profoundly 
affect the asymptotic kinetics. 

To see this explicitly, we solve for the density of dimers 
to find 



C2 = c 2 (0) r e 



2 „-A 
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where we define A = f*dt'a(t'). Using Eq. © we express 
A in terms of r: 



A = A(r) 
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where we further define r = ci(0)/i?(0). 

Depending of the value of 7 we obtain three distinct 
long-time behaviors: 



• For < 7 < 2, A — > Tjz^ — * const as t — > 00. 
Thus the dimer concentration has a universal t~ 2 
asymptotic decay; only the amplitude of the decay 
depends on the initial state. 

• For 7 > 2, the dimer density decays as a stretched 
exponential e~ A times a power-law prefactor. 

• For 7 = 2, the dimer concentration asymptotically 
decays as t~( 2+r ', where the exponent depends on 
the initial concentration ratio r = ci(0)/i?(0). 



III. ASYMPTOTIC MASS DISTRIBUTION 

To determine the asymptotic cluster concentrations C/., 
we re- write the rate equations (JSj in the form 



where 



a 



Ck + ack = at, 



Ci+jR, a k = c 1 c k _ 1 
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Due to the recursive nature of Eq. 
obtain the formal solution for c k : 



Cfe(t) = £(t) \ c fe (0) 



we immediately 

, a k {t')' 
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(10) 



Therefore, 



, Cfc_i, we can compute a k and then 



with £{t) = exp 

once we know c\ , . 
determine c k . 

Before proceeding to general k, it is instructive to com- 
pute C3. From Eq. (| 1 0|) . and using the previous results 
for ci and C2, we obtain: 



c 3 = [c 2 (0)A + c 3 (0)] T- 2 e- x 



(11) 




As in the case of dimers, there are three distinct asymp- 
totic behaviors for the trimer concentration that are best 
expressed in terms of C3/C2. We obtain 



(12) 



that is, C3/C2 oc A for any value of 7 (see Eq. (J7J. While 
the exact expressions for the concentrations c k become 
unwieldy as k grows, the asymptotic expressions turn 
out to be simple: C3/C2 — * A, c 4 /c 2 — > A 2 /2, etc. We 
will confirm this hypothesis by induction in the following 
subsections. 



A. 7 > 2 

When 7 is strictly greater than 2, we easily find that 
Ck/c-2 oc A fc_2 as t — > 00 by explicit calculation for small 
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k. We therefore seek an asymptotic solution for Ck of the 
form 



C2 



k-2 



(13) 



where . . . denotes subdominant terms in the limit t — > oo 
and the amplitude will be determined below. 

For concreteness and simplicity, let us consider the bi- 
disperse initial condition in which only ci(0) and C2(0) 
are non-zero while Cfe(O) = for all fc > 2. To extract 
the asymptotics of from Eq. I|10|) . only the second 
term on the right-hand side is needed. Let us suppose 
that Eq. ljr3*|) holds for C2, . . . , Cfe. To leading order, we 



have otk+i = c\Ck, and substituting Ck 
Eq. i(T0"|) . we obtain 
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into 
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where to obtain the second line we use £{t) 
Thus = l/(fc — 2)! and we arrive at the remarkably 
simple asymptotic solution 



Cfc+2 = c 2 



A^ 



c 2 = c 2 (0) t 2 e 



(15) 



with r and A given by Eqs. Q and Q. Explicitly, we 
have 



Cfc+2 



1 

fc! 



t -2+*(i-a/7) exp ( 



const, x t 



l-2/ 7 \ 
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Thus Cfe has a very different asymptotic behavior than in 
the constant-kernel system 0- As discussed previously, 
this anomaly arises because monomers are very long-lived 
and these then strongly influence the long-time kinetics. 

One important caveat, however, is that the derivation 
of Eq. I |l 5(1 applies only for finite fc. For sufficiently large 
fc the neglected subdominant terms accumulate to pro- 
vide a relevant contribution. The simplest way to see 
that Eq. 1(15(1 cannot hold over entire mass range is to 
use this equation to compute the total cluster density. 
This gives J2k>2 °k — °2 e A = C2(0)r~ 2 , in contrast to 
the correct result i?(0)r _1 given in Eq. Q. 

We therefore conclude that the mass distribution nat- 
urally divides into a small-mass boundary layer, fc < fc*, 
that contains an asymptotically negligible fraction of the 
total mass, and a remaining bulk region. The extent of 
the boundary layer may be determined as the maximum 
of Ck in Eq. and gives fc* = A oc i 1-2 / 7 . Within 
this boundary layer, the mass distribution grows rapidly 
with mass, but this region contains only an asymptoti- 
cally negligible fraction of the total mass. On the other 
hand, most of the mass, which exhibits conventional scal- 
ing behavior is contained in the bulk region. 

In the marginal case of 7 = 2, the asymptotic solu- 
tion of Eq. I(15|) is still valid. Since now A = r lnr, 



the density of dimers decays algebraically in time as 
C2 = 02(0) r~( 2+r ), and the width of the boundary layer 
grows logarithmically fc* = A = rlnr. The following 
main part of the mass distribution is again characterized 
by conventional scaling behavior (see Eq. 1(1911 below). A 
justification of this picture is given in Appendix. Thus 
for 7 > 2, the mass distribution exhibits two growing 
scales and conventional single-mass scaling is violated. 



B. 



7 < 2 



By solving c 2 ,c 3 , . 
esis that for fc > 2 



exactly, we are led to the hypoth- 



ec oc r 



(17) 



This indeed can be directly checked by induction. More- 
over, in the scaling limit, 

fc 
1 



00, t 



— = finite. 



(18) 



the mass distribution admits the conventional scaling 
form 



7 2 £ 2 



-2fc/( 7 i) 



(19) 



Summing ^ Ck we indeed recover R = 2/ (ft), while the 
next moment kck equals one (based on the total mass 
set equal to one initially). Equation ((19(1 also describes 
the scaling portion of the mass distribution when 7 > 2. 

The origin of the conventional scaling behavior for the 
case 7 < 2 is simple. For 7 < 2 monomers disappear 
quickly, since c\ oc £~ 2 / 7 , and their asymptotic influ- 
ence is negligible. Thus the reaction effectively reduces 
to a constant-kernel system that begins with dimers. The 
monomers do influence the small-mass behavior, namely 
Cfe = Akt~ 2 with mass-dependent amplitudes Ak] e.g., 
A 2 = c 2 (0)[7i?(0)/2]- 2 e- A and A 3 = A 2 [A+c 3 (0)/c2(0)]. 
However, as fc grows the amplitude Ak approaches the 
constant value 47~ 2 . 

The existence of scaling can be proved rigorously, e.g., 
by the generating function approach given in the Ap- 
pendix. A simpler approach is to merely assume that 
scaling holds and check its consistency a posteriori. In 
the continuum limit, we have checked the correctness of 
Eq. 1(19(1 . which is the conventional scaling form for con- 
stant kernel aggregation Q- 



IV. REACTIVE MONOMERS 

We now assume that monomer-monomer reactions do 
occur: Ku = e > 0. This situation has already been 
studied by Hendriks and Ernst The rate equations 
for monomers and dimers now read 



ci = —c\{ec\ + R) 

C2 = ^cf - c 2 (ci + jR), 



(20) 
(21) 
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while the rate equations for clusters with k > 2 are the 
same as in Eq. 0. Correspondingly, the density of reac- 
tive clusters R evolves according to 



Note also that (// — 1)£?2 = eB 2 /2. To solve the recur- 
sion (|26[l we introduce the generating function B{z) — 
J2k>2 Bkz k to reduce to 



R 



- 7 R 2 



(22) 



It does not seem possible to solve the coupled non- 
linear equations (|20|l and l|22|) for ci and R explicitly. If 
we treat C\ as a function of R, however, we can reduce 
Eq. H20|) and i|22[l to a single differential equation that 
can be solved to establish a functional relation between 
ci and R Q]. From this relation, we can deduce that 
different behaviors emerge depending on whether 7 < 2 
or 7 > 2. 

A simpler and more fruitful way to proceed is to seek 
the asymptotic behavior of ci and R without explicitly 
solving the governing equations. We anticipate that there 
are three possible asymptotic behaviors: (i) c\ <C R, (ii) 
ci oc R, and (iii) ci 3> R. Substituting each of these 
scenarios into Eqs. (|2U[1 and (12211 shows that the latter 
case is impossible. Let us now analyze the first two pos- 
sibilities. When ci < R, Eq. (JUJ gives R ~ 2/( 7 f), 
and then Eq. <|2T?|l yields Ci oc i -2 / 7 . Thus the relation 
ci <C R holds when 7 < 2. In the complementary regime 
of 7 > 2, we find that c\ oc R is consistent; furthermore, 
Eqs. (|2njl and (O now give R oc c\ oc t~ x . 

We now analyze these two cases in more detail. 



A. 7 > 2 

When 7 is strictly greater than 2, we substitute the 
ansatz 



(23) 



into Eqs. I|20|l and l|22|) to get a quadratic equation for 
the amplitude A. The physical requirement that A and 
B are both positive fixes the solution to be 



A _ e-l + y /e(e + 1 -2) g = 1-A 
57 — 1 e 



(24) 



The singularity at e = 7 _1 is only apparent and may be 
resolved by applying the 1'HospitaPs rule, to give A — 
7/[2( 7 — 1)] in this case. 

From the formal solution we can check by induc- 
tion that for finite fc, 



c fe ~ B k r 



(25) 



The amplitude B k is found by substituting the ansatz 
(|2*5|l into to give the recursion 



(ji - l)B k = BB k -i + I B * B 

i+j=k 

for k > 2, where we define 

H = j A + B = 1 + v /l + ( 7 -2)/e 



(26) 



(27) 
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Bz 



H-l-Bz 



When k ^S> 1 but still within the boundary layer, the 
asymptotic behavior of B(z) leads to the amplitude 



B k ~ Ck~ 3/2 f3 k 



with 



/3 = B 



C = 



/ey 



7V27T V 1 + y/ey 



(28) 



(29) 



One can verify that j3 < 1 for all 7 > 2 and e > 0; there- 
fore, apart from a power-law prefactor, the amplitude B k 
decreases exponentially with k. 

As in the previous section, the results l|2*5|) and f2"%|l 
are valid only within the boundary layer. Thus the total 
density of clusters within the boundary layer is asymptot- 

J2 k>2 B k =B(z = l), 



ically equal to A't 1 , where A' 
or 



A 



, _ V - 1 - B - y/Qi - 1 - B) 2 - ejB 2 
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(30) 



Since the total density of clusters scales as At 1 , only a 
fraction A' /A of clusters lies within the boundary layer. 
Using Eqs. and CP, we find A - A' = 2/7; this 
implies that the density of clusters in the bulk decays as 
2/(7^). The same holds for 7 < 2, although in this case 
the density of clusters in the boundary layer is asymp- 
totically negligible. 

To determine the asymptotic behavior in the bulk (k — > 
00 and t — > 00), we simplify the rate equations by neglect- 
ing the subdominant terms c\(c k ^x — c k ) on the right- 
hand side of Eq. 0. Further, the sum in has three 
contributions. When i lies within the boundary layer, we 
can replace £\ CiC k -i by c k ^ B^ 1 = CkA't^ 1 ; a sim- 
ilar contribution arises when j lies within the boundary 
layer. In the bulk we replace the sum by the integral. 
Thus 



t+j=k 



jA' 



■ Cfe . 



Combining all terms and using the identity j(A—A') = 2 
we finally convert Eq. to 



' Cfe 



di CiC k - 



(31) 



This equation appears in the standard constant-kernel 
aggregation and its solution, satisfying the aforemen- 
tioned conservation laws J dk c k = 2/7 and J dk kc k = 1 , 
is given by Eq. (|19|) . 
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By matching the mass distribution in the boundary 
layer (c^ oc t^ 1 /3 fe ) and in the bulk (c^ oc t~ 2 ), we esti- 
mate the width of the boundary layer as 



k* 



hit 



ln(l//3) 



(32) 



with P given by Eq. (|29|) . This is a much slower growth 
than in the case of sterile monomers where k* oc t 1 ^ 2 /"' . 

When 7 = 2, we again expect marginal behavior. Pro- 
ceeding as in the case of 7 > 2, we find 



Cl 



(33) 



From the formal solution (|10J) we can check by induction 
that, 



c k ~B k t- l {t\nty 



(34) 



For k > 2, the amplitude B k is found from the recursion 
ffity with = 2, B = 1, and B 2 = e/2. The generat- 
ing function and the coefficients -B^ and C are given by 
taking the corresponding formulae in the 7 > 2 case and 
setting 7 = 2. 

This behavior holds in a boundary layer whose width 
now grows as 



k* 



hit 



ln(Alni) 



A 



1 + V21' 



(35) 



while in the bulk we recover the ordinary scaling mass 
distribution c k = t~ 2 e~ k l l . 



B. 7 < 2 

We already argued that R ~ 2/(-yt) and c\ oc r 2 /^ for 
7 < 2. Solving then for dimers we obtain: 



c 2 oc 



f 



t l-4/ 7 

Generally, we find 

r t- 2 



Cfc oc < 



t~ 2 lnt 

.fe-l— 2fc/7 



< 7 < 4/3, 
4/3 < 7 < 2. 



0<7<^- 



2fc 

fc+1 • 



(36) 



(37) 



^T<7<2. 



Thus for every 7 < 2, small-mass clusters have abnor- 
mal kinetics: only monomers for 7 < 4/3, monomers 
and dimers for 4/3 < 7 < 3/2, monomers, dimers, and 
trimers for 3/2 < 7 < 8/5, etc. The remaining clus- 
ters decay as t~ 2 , and conventional scaling describes the 
cluster-mass distribution. 



V. SUMMARY AND DISCUSSION 

For irreversible aggregation with distinct monomer- 
monomer (Kmm), monomer-cluster (.Kmc); and cluster- 
cluster (Kcc) reaction rates, the dynamics depends cru- 
cially on the ratio 7 = Kqc/Kmc, while e = ^mm/^mc 
plays a lesser role — all that matters is whether e = 
or e > 0. For e = and 7 < 2, there is conventional 
scaling with a single mass scale growing linearly with 
time. For 7 > 2, there are two scales: the boundary layer 
k < k* oc £ 1-2 / 7 where the mass distribution has an un- 
usual Poisson form, and the bulk region k > k* where 
conventional scaling holds. When e > and 7 < 2, there 
is conventional scaling, except that light clusters have 
abnormal kinetics - monomers for 7 < 4/3; monomers 
and dimers for 4/3 < 7 < 3/2; monomers, dimers, and 
trimers for 3/2 < 7 < 8/5, etc. When e > and 7 > 2, 
the behavior in the boundary layer is very different from 
that when e = 0. In particular, the mass distribution 
decays with mass while for e = the mass distribution is 
peaked. 

A possible application of this model is to clustering on 
surfaces. Consider a two-dimensional substrate with dif- 
fusing single-layer islands that aggregate whenever they 
meet. In the diffusion-controlled limit, the reaction rate 
of an island of radius R and diffusivity D is proportional 
to D In R [j| . To a good approximation, we can ignore the 
island radius and think of point-like islands that always 
occupy a single lattice site. When such a cluster hops 
onto already occupied site, two clusters immediately co- 
alesce into a single cluster. 

If monomers hop with rate D while all heavier clusters 
hop with the same unit rate, then the corresponding re- 
action rates for point-like islands are K\\ = 2D, K\j = 
1+D, Ka = 2. As a result 



K11 
Id, 



2D 
1 + D' 



K 



1 + D 



(38) 



For immobile monomers, D = 0, we have e = and 
7 = 2. In this case, the cluster mass distribution is 



c k+ ,=c 2 (Q)r- {2+r) \r 
k\ 



(39) 



within a boundary layer that grows logarithmically with 
time k* = rln(l + R[Q)t), where r = ci(0)/iZ(0). In the 
bulk, ordinary scaling holds in which c k = t~ 2 e~~ k l l . 
When D > 0, 7 is always less than 2. Therefore, 



R = 



l + D 



t 



Ci oc t 



-1-D 



(40) 



The decay rate for clusters of mass k > 1 is also simple: 



~ 2 D>l/k, 

c k oc { r 2 hit D = l/k, 
t -i-kD D < 1 j k _ 



(41) 
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APPENDIX A: GENERATING FUNCTION 
APPROACH 



The term t" 1-1 ^ 1- ^ in the denominator can be ig- 
nored within the boundary layer. Therefore, the denom- 
inator becomes (1 — z)[l + (l + r) _1 z], and the generating 
function simplifies to 



G(z,t) = c 2 {t)z 



2 \z 



2 + r 



l + r 



l+r 



(A6) 



In the case of sterile monomers, we can employ the 
generating function approach as an alternative to a direct 
solution. The generating function 



G{z,t) = Y,z k c k {t) 



(Al) 



k>2 



recasts an infinite set of rate equations into a single dif- 
ferential equation 



dG 7 2 

~dt ~ 2 +(^-«) G - 



(A2) 



This is a Bernoulli equation that is readily solved in terms 
of G -1 . The solution is 



G = 



r 2 p (l-*)A 



T dr' 



'LV e (i-*)(A-V) 

G(z,0) J l R(0) \t'J 



(A3) 



where A = A(r) and A' = A(r'). For the bi-disperse initial 
condition, Cfe(O) = c\ (0)<5fc,i+C2 (0)5k,2 : we have G(z, 0) = 
C2(0) z 2 and R(0) = c 2 (0). Using these relations together 
with C2 = C2(0)r _2 e _A we re- write l|A3(l as 



er zX e- zX f T dr' 



c 2 z z 



C-2 



-(l-z)A' 



1 {r'f 



(A4) 



We now consider in detail the marginal case of 7 = 2 
and justify the two-scale structure of the mass distribu- 
tion. When 7 = 2, we have A = rlnr and (|A4|> becomes 



G(z,t) = 



c 2 (t) z 2 e Xz ( 



1- IT"* 
l+r 



l+r l+r L 



(A5) 



To extract the mass distribution, we expand G(z,t) in a 
Taylor series in z. This gives 



Cfe+2 
C2 



k 

E 

n=0 



A™ 
n! 



1 



1 



1 



■ (A7) 

2 + r 2 + r V l + r/ 
Thus apart from the leading contribution that equals 
A fe /fc!, in agreement with our previous result in Eq. Ijl5|l . 
we find all the correction terms; e.g., the leading correc- 
tion is (1 + r)- x \ k - 2 /{k - 2)1 

When both the quadratic and the transcendental terms 
in the denominator balance each other, the generating 
function accounts for the main part of the mass distribu- 
tion. The both terms are comparable when 1 — z oc t . 
Writing 1 — z — £/t and taking z — > 1, r — > 00 limit, 
with £ kept finite, we simplify the generating function: 



G 



c 2 (0) 1 
t l + (2 + r)C 



(A8) 



This form of the generating function implies that 



c 2 (0) 
(2 + r) 7 



exp 



1 k 



2 + r t 



(A9) 



Since 7 = 2 and (2 + r)R(0) = 2c 2 (0) + ci(0) = 1 (for the 
bi-disperse initial conditions i?(0) = 02(0) and the mass 
density is always set equal to one) , the above scaling form 
reduces to the anticipated scaling form of Eq. (tT^ll . 
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